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Instability of 6 dimensional Rubakov-Shaposhnikov model is reinvestigated. It is shown that the 
model is unstable in scalar perturbations sector with very particular instability pattern: there are 
no unstable modes for the first two lowest angular harmonics, m = and m = 1, whereas there is a 
single negative mode for each higher m. 
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I. INTRODUCTION 



In the beginning of 80th two very important papers by Rubakov and Shaposhnikov were published in the same 
issue of Physics Letter B [l[ , . The first one p| was discussing possibility that (in modern language) we live on a 
brane in higher dimensional space and in the second one Q the warped compactification was introduced in order to 
attack cosmological constant problem. These ideas created basis for "extra dimensional revolution" which happen 15 
years later @, gj, 1, 0,0,1- 

The stability of the Rubakov-Shaposhnikov model with warped compactification Q was questioned [0] soon after 
the model was suggested. It was found that the model is stable under tensor and vector perturbations, but has 
unstable modes in scalar perturbations sector. Recently we became aware [l4| that there is an algebraic error in the 
prove namely in the Eq.(23) of @. The aim of present note is to to correct this error and reinvestigate the stability 
of the Rubakov-Shaposhnikov model. Since in addition there are numerous misprints in the most of equations in the 
journal version of [9j, first we repeat here derivation of Schroedinger equations and then give direct numerical proof 
of existence of an unstable modes in scalar perturbations sector. 



II. SPONTANEOUS COMPACTIFICATION WITH ZERO COSMOLOGICAL CONSTANT 

The solution, leading to zero four dimensional cosmological constant, proposed in Q looks as follows. We consider 
gravity in d + N— dimensional space-time with the metric ()ab (signature H — ... — ). The Einstein equations is written 
with the cosmological constant: 

Rab ~ ^9abR = ^9ab ■ (1) 



It is assumed that A > 0. With the warped ansatz for the metric 

g ab (x a ) 



SAB= (*(*)9 r V) J , ( (2) 



the Eq. JTJ) reduces to the system of equations 



R/MV t^Qla/R — Aphysff/ii; j (3) 

R ab = Ag ab + d(M^ _ gsggg) , (4 ) 

ab N + d~2 y 2a 4cr 2 ' K ' 

iv.V'a + (d - 2)X^ — — o~ = . (5) 

2 V ; 4ct N+d-2 d-2 y ' 
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the hats and tildas respectively denote d + TV-dimensional and iV-dimensional quantities, p, v, ... = 0, 1, d — 1, and 
a, b,... = d,...,d + JV- 1. 

In these equations A p h ys is an arbitrary parameter arising from the separation of variables, is constructed from 
g^v and and R au and V a are constructed from g^ according to the usual rules. The Latin indices a, b, ... are raised 
and lowered with the metric g au . Note that, in contrast to the standard approach to spontaneous compactification, 
the space defined by the metric c/ab is not the direct product of the d- and TV-dimensional spaces. This difference is 
related to the presence of a warped factor a{x a ) in front of g^, v {x) and is critical future in the entire discussion. 

The Eq. ([3]) is the Einstein equation for the c£-dimensional metric g^ with the cosmological constant A p h ys . For 
Aphys = it has a solution corresponding to a flat space. Assuming d — 4 and N = 2, the equations (jlj and (0 can 
be solved with the result: 



■Jul, 



-1 
f(p) 



f(p) 
<r{p) 



5A 



,5A , 
tg(\/-5-P) 



2 r 



/5A 

cos(y — p) 



-i 4/5 



5A , 
cosh/ —p) 



4/5 



(6) 
(7) 

(8) 



where A ph ys = 0, x 4 = p, p € [0, p max ], x 5 = 9,0 E [0, 27r]. 

In spite the fact that this solution is noncompact in the usual sense (the circumference x^ — const, p — const can 
be arbitrarily large, I = 27TV— J — > oo as p —> p ma , K ), it can be shown that the presence of the two extra dimensions 
is unobservable at low energies 0. 



III. STABILITY ANALYSIS 



A. The equations of motion 



Let us find the equations of motion for a small fluctuations about the solution g AB . For this we substitute (jab 
9ab + ( ab into the Einstein equations Eq. ([I} . Neglecting all terms with powers higher than first in cab , we find 



VcVacb + VcVse^ - V c V c e AB - V ' A V Bee = 



4A 



2-N-d 



CAB 



(9) 



where the covariant derivatives Va are calculated using the background metric (Jab- l n wna t follows we set d = 4 and 
N = 2. Next we use the fact that g AB is independent of x^ and go to the momentum representation in a;„. Denoting 



and substituting the decomposition of /i M „ and A afl into components with spin 0, 1 and 2 



(10) 



(11) 



where 



#V=0, fc^r aM = 0, fc^ = 0, I/ = 0, (12) 
from the Eq.Q we obtain seven separate equations for tensor, vector and scalar (under rotations of x^) perturbations: 

k 2 a~\ u - (V a V a + 27r a V a )^ = ; (13) 



i(V a + n a )r av - a(V a V Q + 2^ a V Q )/„ = , (14a) 
k 2 a- x r all + ik 2 \7 a ff, - (V h V b + 7r b V h - A)r aM + (V 6 V a - 7r a V fc + 27r fc V a - V b 7r a - 2n a TT b )r^ b = ; (14b) 
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k 2 a- x S + ik 2 a-\ a $ a ~ ^7r a V a (P + 35 + tp b b ) + (2n a n b + VV + n b V a )tp ab - (V a V a + 27r a V a )5 = , (15a) 

2ik 2 a~ 1 {V a + Tr a )<t>a + k 2 a-\2S + <p b b ) - (V a V a + 27r a V a )(P - 5) = , (15b) 
i(V b ^a b ~ V a (35 + v b b ) + \naVb + *b<p a h ) - (V b V b + 2v b V b - A)$ a 

(VfcVa - TTaVft + 27T fc V Q - VfeTTa - 27r a 7r 6 )$ h = , (15c) 
feV^Vab + ifc a £T _1 (V tt *6 + V 6 $ a ) - -(WftVa + 7T a Vfc)(P + 35) - V a V b (P + 35) 

V a V^ c c + (V c + 27r c )(V a v?f, c + V b ^ a c ) + (A - V C V C - 2TT c \7 c )ip ab = , (15d) 



where 7r a = (d/dx a )\na. 

The system of equations f|13til5[) is a system of eigenvalue equations with the role of the unknown eigenvalue played 
by k 2 . If a system (for example, (|15p) is consistent, there will be at least one equation of the form 

LV(k,x a ) = kH(k,x a ) , (16) 

for some ^(k,x a ), where L is a differential operator acting on x a . The question of stability of the solution g° AB now 
reduces to the question of the possible values of k 2 . The solution will be linearly stable if k 2 has no negative values 
and solution is unstable if k 2 can take negative values. 

Note that in this problem there is gauge invariance related to the invariance of the equation ([1]) under general 
coordinate transformations x A — > x' = x A — Ax A . In terms of variables Eq. (|10[) and Eq. dill) , gauge transformations 
with the parameters Ax A = [o p_1 (a; A ' + k^r]),x a ] looks like 

d d 

6ln„ = Q, 5r a ^, = (^-^ - 7r )a; M , 6f ll = ia~ u^, 6<$> a = - 7r a )?7 + ixa, 

5S = TT aX a , 6P^2ik 2 <J- 1 7 1 + TT aX a , Sl Pab = + gac-^ + Qbc^X* ■ (17) 

B. The mass spectrum 

It follows directly from the Eq. ([13]) for tensor perturbations that k 2 is non-negative for this sector. In order to see 
this, we set fi, v = in this equation, multiply both sides by a 2 \/—fl 00 (—k) — [— det^^g)] 1 ^ 2 ! 00 ^ k), and integrate 
over dx^dx 5 = dx. Integrating the right hand side by parts (the correctness of this procedure can be rigorously 
justified), we obtain the equation 

a ;00/ 



aV-fdx\l 00 (W = -J <TW-fdx[V a l 00 (k)}* [V a / uu (fc)] , (18) 

from which it follows that k 2 > 0. 

For vector perturbations first we have to fix the gauge. We choose the gauge conditions as 

U=0. (19) 

Furthermore, we set /j, = in the second of equations (|14p. multiply by a\/—fr a0 (—k), sum over a, and integrate over 
dx. Then with the help of gauge condition Eq. (flU]) and use of the first of equation (TTi|) after some transformations 
we find 

fdxr a0 r aQ = J a^f~fdx{-V a r m V a r w - \ a i: b r aa r b - ^K c n c r a \ a0 ) . (20) 

The integral multiplying k 2 and the right hand side are both non-positive. Therefore, for the vector perturbations 
also k 2 > 0. 

Let us now turn to equations Eq. (|15p for scalar perturbations. We choose the gauge condition in the form 

$ a = 0, P-S = 0. (21) 

We shall assume that k 2 ^ 0. Then equation (fl5k) is a consequence of the three other equations and can be omitted. 
From Eq. (|15b ) we obtain 

S=~\ipb, (22) 
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which can be taken as a definition of S in terms of ip a b - Using this equation and the gauge condition Eq. (j2"Tj) . the 
two remaining equations can be written in terms of <p a b '- 

(V a + Tr a )ip b b = -2(V 6 + 7r b )(^ a b , (23a) 
k 2 a- 1 ip ab = -(V a V b + 7r Q V h + ir b \7 a )tp c c - (V c + 27r c )(V a ^ c + V fcV 3 a c ) + (V C V C + 2tt c V c - A)cp ab . (23b) 

Equation fZ3k ) gives two relations between three variables <f> ab , so a single independent variable remains. Three 
equations (|23b ) are equivalent to each other and determine the spectrum of k 2 . The problem is to solve the constraint 
(|23k ). that is, to express all the in terms of a single independent variable £ and its derivatives. Then from (|23b) we 
find following equation for £: 



fc 2 £ = M£ , (24) 



where M is a differential operator acting on x a . 
Let us expand tp a b(k, p, 0) in a Fourier series in 



Vab( P ,0)= J2 ^Vab™^), (25) 



m— — oo 



and consider cases m = and to =^ separately. 
For to = the Eq. (|2"3"r ') gives: 

(9p+r + 7r 4 )( y 5 5 4 = 0, (26a) 
(0„ + ^ 4 )(3^4 4 + V5 5 ) = -2r(^ 4 4 - <^ 5 5 ) , (26b) 

where T = r| 5 = i<9 p ln/. The first of these equations has no regular solutions except for zero, so it gives tp b A = 0. 
Using this condition, the Eq. (|26b ). the background equation ([5]) and denoting 3<y9 4 4 + ^5 5 = £bj we find 

fcV-^o = -a p 2 Co + (2^ - r - 7r 4 )d p £ + (2tt 4 ^ - 2^ 4 r - 7r 4 2 )£ • (27) 
Multiplying this equation by £0 and integrating over p from zero to p max with the weight W(p) given by the expression 

W(p) = exp (- J dp(2^- - r - tt 4 )) > , (28) 



after integration by parts we obtain 



A; 



J o-^W^dp = J Wdp 



(d p £o) 2 + (2tt 4 ^ - 27r 4 r - 7r 4 2 )£ 2 



(29) 



The positivity of the last term in the integrand follows from the explicit form of background solution, (J7]) and (J5J. 
Therefore, k 2 > for to = 0. 

In the to 7^ case we denote (p^ m {p) — ■^■y 5 4 „ l (p). Dropping the index to on the fields f a b: we rewrite Eq. (|23a ) 

as 

(0 P + ^ 4 )(3<^ 4 4 + <^ 5 5 ) = 2^ 5 4 - 2r(< /34 4 - <^ 5 5 ) , (30a) 
<p 4 4 + 3(^ 5 5 = -2(a p + r + 7r 4 )(^ 5 4 , (30b) 
where /1 = —to 2 //. The solution of this system has the form: 

£ = 3^ 4 4 + <p 5 5 + 2r^ 5 4 , G = /i + d p T , 

, 4 (gg + r + 7r 4 )e 4 5 T . (a P + r + 7r 4 )e 

^ - ^ = ^(2m^ 5 4 - (d P + ^ 4 )(3^ 4 4 + ^ 5 5 )) . 
After some awkward algebra, the Eq. (|23b ) gives following for £: 

^-^ = -^ + ^(^ + 5^, (32) 

A( P ) = 2^-r, B(/9) = /i + I^ + 27r4r + 2(r + 7r4) ^ (33) 
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It is convenient to change to the dimensionless variable x — y ^fp, x £ [0,7r/2], in the Eqs. f[32| [33]). 

Now we bring Schrodinger Eq. (|32p into standard form in two steps. First we change independent variable from x 
to r according to dr = dx/y/a. When x € [0, 7r/2] variable r is changing from to T max , 

dx 7T 3 / 2 CS c(^) , . 

Tmax = / 7 ; — TToTF = — ; 7 ; a, ~ 2.27221542 . (34 

J [co8(x)]*/* 2r(X)r(|) v ' 

This way we get rid of factor a in the l.h.s. of Eq. (|32|) and coefficients ^4 and B are changed to 



A r = y/aA + ^y/an 4 ,B T = <tB . (35) 
Next with the transformation £ = exp(i J A T dr)x we get rid of the first derivative term and arrive at 

(—^ + U)x = k 2 x , with U = B T + ^A 2 T -^d T A T . (36) 



Potential U close to r — > behaves as 

4 T 2 



35 

U = — + 0(t 2 ) for m = 1 , (37) 



and 



and close to t„ 



i 



56 — 8m 



U= ' 4 + + 0(r 2 ) for m>2, (38) 

r z 15 



U = - Mt _ r)2 - 09) 

Note that quantum mechanical potential U = —j/x 2 for 7 > 7 cr = 1/4 corresponds to unstable situation ("falling" 
to the center, see e.g. [HI)- So, in Eq. (|36l l39| exactly critical case is realized, which is on a border between stability 
and instability. Another observation is that the m = 1 case is distinguished, because potential U is not negative in 
the inner region, while starting from m = 2 it is negative not only asymptotically r — > r max , Eq. (|39[) . but also in the 
inner region, Fig. 1. 

The regular branch of the wave function \ close to r — > behaves as 

x oc t 7 ' 2 - -V 11 / 2 for m = l, (40) 



and 



and close to r — > r mnT as 



16 



xocrm +i/ 2 + 56 8m 2 15fc 2 Tm+5/2 for to ^ 2j 



60(m+l) 



X^lw-r) 1 / 2 . (42) 

To determine the number of bound states of Schrodinger equation in a given potential we investigated the zero energy 
wave function. According to known theorems (see e.g. fllj ]) the number of nodes of zero energy wave function exactly 
counts the number of negative energy states. Solving numerically Schrodinger Eq. (|36p with above boundary conditions 
we found that there are no negative modes in m = 1 case, whereas there is single negative mode for each higher m. 
We checked this statement up to m — 10. 



IV. CONCLUDING REMARKS 



We have shown that the solution described by Eqs. (|2I6I7I8P and corresponding to the A p h ys = is linearly unstable. 
The instability is related to ^-dependant perturbations, which are scalars under rotation of the four dimensional 
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coordinates x^. We found that a single unstable mode appears in spectrum of linear perturbations for each angular 
harmonic with m > 2. A similar situation can be expected to arise for small A p h ys - Even if the solution is stable 
starting from some A° h ^ 0, it is quite improbable that the value of Ap hys will be ~ 10 _56 cm~ 2 , in agreement with 
current observations [12|, [l3j . 

Although the Rubakov-Shaposhnikov solution is found to be classically unstable, knowledge of unstable modes 
can be useful, since it suggests the form of the stable solution to be sought. Since the perturbations leading to 
instability are asymmetric under ^-rotations, it is clear that the initially symmetric state of the system tends to the 
more favored asymmetric state. Lorentz invariance is not violated in the development of the instability, as it would 
be if the instability were related to the vector perturbations. So, it is quite possible that the Eqs. (14151) for N = 2 
have asymmetric, ^-dependent solution, which might be stable. 

Since Rubakov-Shaposhnikov model is basic ingredient for many modern higher dimensional setups it is natural to 
ask whether the instability disappears by adding extra fields. So, question of stability should be carefully checked in 
each case. 
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FIG. 1: Shape of the potential U(t) for m = 1, 2 and 5. 



